Generalized Cloaking and Optical Polyjuice 
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In 2006, it was suggested by Pendry et al. [1] that an object surrounded by a coat- 
ing consisting of an exotic material becomes invisible to electromagnetic waves. This 
device was named "invisibility cloak" in reference to Harry Potter, the popular char- 
acter of J.K. Rowling. Beside his famous cloak, the little wizard has other spells to 
go unnoticed. Among the most spectacular is the "polyjuice potion" that is able to 
turn somebody into anybody else's appearance [2]. In this paper, we do not present 



a potion but rather an optical device able to accomplish the same task, i.e. to give 
an arbitrary optical response chosen in advance to any other object placed inside the 
device. In fact, the principle is here very similar to the design of Pendry's invisibility 
cloak but, instead of geometrically transforming an empty domain, it is a region con- 
taining the object to be imitated that is transformed leading thus to a generalization 
of cloaking. 

In this section, we present a generalization of cloaking able to arbitrarily transform 
the electromagnetic appearance of an object. The basic principle is to obtain the 
constitutive relations of the cloak by application of a space transformation to a non- 
empty region. 

A geometric transformation is given by a map <p from a space N to a space M. For 
all our practical purposes, M and N will be here the whole or parts of M 3 . Given a 
Cartesian coordinate system x on M and an arbitrary coordinate system x' on N, (p : 
N — > M is described by x(x'), i.e. x given as function of x'. All the useful information, 
i.e. necessary to transform differential forms (indeed the formalism of differential 
geometry is the most natural to write the Maxwell's equation [3-6] specially when 
quite general geometric transformations are involved) and other covariant tensors, is 
contained in the Jacobian matrix field J(x') = <9x(x')/<9x'. 

The basic principle of transformation optics is that, when you have an electromag- 
netic system described by the tensor fields e(x) for the dielectric permittivity and 
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/x(x) for the magnetic permeability in the space M, if you replace your initial mate- 
rial properties by equivalent material properties given by the following rule [7-10]: 

^(x') = J- 1 (x')|(x(x , ))J" T (x , )det(J(x / )), 

^(x') = J- 1 (x')^(x(x'))J- T (xO det(J(x')), (1) 

(J~ T is the transposed of the inverse of J), you get an equivalent problem on N. 
Here, an equivalent problem means that the solution of the new problem on N, i.e. 
electromagnetic quantities described as differential forms, are the pulled back of the 
solution [5] of the original problem on M and that the same Maxwell's equations 
(i.e. as if we were in Cartesian coordinates or, more accurately, having the same form 
written with the exterior derivative) are still satisfied. 

In the case of the cylindrical Pendry's map [1,11,12], described by the transforma- 
tion of the 2D cross section, the plane IR 2 minus a disk Di of radius R\ is mapped 
on the whole plane M 2 in such a way that a disk D 2 of radius R 2 > R±, concentric 
with Di, is the image of the annulus D 2 \D\ by a radial transformation. In cylindrical 
coordinates, this transformation is given by: 

r = (r' - R 1 )R 2 /(R 2 - R x ) for R l < r' < R 2) 

(2) 

= ff, z = z'. 

As for the outside of the disk D 2) the map between the two copies of 1R 2 \-D 2 is the 
identity map. 
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The material properties given by rule ([I]) corresponding to this transformation 
provide an ideal invisibility cloak: outside D2, everything behaves as if we were in 
free space, including the propagation of electromagnetic waves across the cloak, and 
is completely independent of the content of D\. 

Now, rule ([!]) may be applied to D2 containing objects with arbitrary electromag- 
netic properties so that a region cloaked by this device is still completely hidden but 
has the appearance of the objects originally in D2. We may call this optical effect 
masking [13] or "polyjuice" effect. 

Figs. [2] and [3] show the effect of masking on a scattering structure. On Fig. [2j 
a cylindrical TM wave emitted by a circular cylindrical antenna is scattered by a 
conducting triangular cylinder (the longest side of the cross section is 1.62A and 
e r = 1 + 4(H). The field map represents the longitudinal electric field E z (x,y) and 
the outer boundary of the cloak is shown to ease the comparison with the masked 
case. On Fig. |3j the same cylindrical TM wave is scattered by a masked triangular 
cylinder (but the diffracting object inside the cloak may be arbitrary as far as it is 
small enough to fit inside the cloak). This triangular cylinder is the symmetric of the 
previous one with respect to the horizontal plane containing the central fibre of the 
cylindrical antenna. This bare scatterer would therefore give the Fig. [2] image inverted 
upside-down but, here, this object is surrounded by a cloak in order to give the very 
same scattering as before. Indeed, on both sides, the electric fields outside the cloak 
limit are alike. 
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Fig. [4] highlights the different scattering patterns by displaying the value of 3te(E z ) 
on a circle of radius 4A located around the antenna-scatterer system in the three 
following cases: the case of Fig. [2] (original) with the triangle alone, the case of Fig. [3] 
(coated), and the triangle of Fig. [3] without the coating (reversed). It is obvious that 
the coating restores the field distribution independently of the object present in the 
central hole. 

The numerical computation is performed using the finite element method (via the 
free GetDP [14] and Gmsh [15] software tools). The mesh is made of 148,000 second 
order triangles including the Perfectly Matched Layers used to truncate the compu- 
tation domain . The singularity of e and \i requires a very fine mesh in the vicinity of 
the inner boundary of the cloak (see Fig. [I]) and is also responsible for the small dis- 
crepancies between the numerical model and a perfect cloak (see Fig. [3]) — including 
the non zero field in the hole of the cloak. 

Note that a small technical problem arises in practice when rule ([!]) is applied: 
the material properties are defined piecewise on various domains and it is very useful 
to know explicitly the boundaries of these domains, e.g. to build the finite element 
mesh (see Fig. [I]). These boundaries are curves in the cross section and are thus 
contravariant objects. Therefore, their transformation requires the inverse map y?" 1 
from M to N. Fortunately, map ^ is very simple to invert. 

More explicitly, for a given curve x(i) of parameter t in the initial Cartesian coor- 
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dinates, its push forward by Pendry's map is: 



x'(t) = ^ 1 (x(t)) = ( 



i?2 — Ri 

R~2 



+ 



x(t)|| 



R 



)x(t) 



(3) 



with the same variation of the parameter t. Note that the most common curves used 
in the design of devices, i.e. line segments and arc of circles, are transformed to less 
usual curves except for radial segments (with respect to the center of the cloak) and 
arc of circles concentric with the cloak. 

On Fig. [3j the image by y? -1 of the triangle of Fig. [2] is the curvilinear triangle 
inside the coating region of the cloak. In practice, this anamorphosis of the triangle 
is described by three splines interpolating each 40 points that are images of points of 
the segments by if' 1 . 

Transformation optics do not offer only the possibility make optically disappear 
objects in invisibility cloaks but also to completely tune their optical signature i.e. 
to give them an arbitrary appearance. The possibility to place an object inside the 
coating of the cloak and that it will therefore appear different was already considered 
in [11] where a point source was shifted creating a mirage effect. A more general 
case is considered here since both the shape and the position of the object placed 
in the coating are modified . Note that if the object used to create the illusion is 
perfectly conducting and surrounds the central point of the cloak, its anamorphosis 
will surround the inner boundary of the cloak and it de facto suppresses the singular 
behavior of the material properties. In this particular case, the present approach is 
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similar to the the hiding under the carpet idea of Li and Pendry [16] but here a 
bounded object replaces the infinite reflecting plane. This technique can be naturally 
extended to cloaks of arbitrary shapes [17]. 
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List of Figures 

1 This figure shows a part of the triangular mesh used for the finite 
element modeling of the scattering problem of Fig. [3j The singular 
behavior of the permittivity and of the permeability requires a very 
fine mesh along the inner boundary of the cloak in order to achieve a 
satisfactory accuracy with the numerical model [11] 

2 A conducting triangular cylinder is scattering cylindrical waves. . . . [£2] 

3 A triangular cylinder different from the one on Fig. [2] is surrounded 
by a cloak designed to reproduce the scattering pattern of the Fig. 
[2] triangular cylinder in spite of the change of diffracting object. Of 
course, the diffracting object inside the cloak may be arbitrary as far 

as it is small enough to fit inside the cloak [T3] 

4 The value of the electric field (the real part of E z ) on a circle of radius 
4A concentric with the cloak is represented as a function of the position 
angle 6 (increasing counterclockwise and with 8 = corresponding to 
the point the most on the right). The three configurations considered 
here are the ones of Fig. [3] (coated), Fig. [2] (original), and the triangle 

of Fig. [3] without the coating (reversed) d3] 
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Fig. 1. This figure shows a part of the triangular mesh used for the finite 
element modeling of the scattering problem of Fig. [3] The singular behavior 
of the permittivity and of the permeability requires a very fine mesh along the 
inner boundary of the cloak in order to achieve a satisfactory accuracy with 
the numerical model. 
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Fig. 2. A conducting triangular cylinder is scattering cylindrical waves. 
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Fig. 3. A triangular cylinder different from the one on Fig. [2] is surrounded by 
a cloak designed to reproduce the scattering pattern of the Fig. [2] triangular 
cylinder in spite of the change of diffracting object. Of course, the diffracting 
object inside the cloak may be arbitrary as far as it is small enough to fit inside 
the cloak 
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Fig. 4. The value of the electric field (the real part of E z ) on a circle of radius 
4A concentric with the cloak is represented as a function of the position angle 
9 (increasing counterclockwise and with 6 = corresponding to the point the 
most on the right). The three configurations considered here are the ones of 
Fig. [3] (coated), Fig. [2] (original), and the triangle of Fig. [3] without the coating 
(reversed) . 
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